AD-A174  265  RELAXATION  FUNCTIONS  MENORV  FUNCTIONS  RANDOM  FORCES  AND  1/1 
ERGOD IC I  TV  IN  THE  <U>  GEORGIA  UNIV  ATHENS  DEPT  OF 
PHVSICS  J  FLORENCIO  ET  AL  05  NOV  86  TR-22 
UNCLASSIFIED  N00014-84-K-0265  F/G  12/1  NL 


flic  FILE  COO  AD-A174  365 


$ 


_ Unclassified _ . 

SECURITY  CLASSIFICATION  of  THIS  FAC®  r*h«i  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


miT-i-i  it 


I.  REPORT  NUMBER 

Technical  Report  No.  32 


4.  TITLE  (md  Subtltlt) 

Relaxation  Functions,  Memory  Functions,  Random 
Forces,  and  Ergodicity  in  the  One-Dimensional 
Spin-li  XY  and  Transverse  Ising  Models 


7.  AUTHORf»J 

Joao  Florencio  Jr,  and  M.  Howard  Lee 


’REFORMING  OROANIZATION  NAME  ANO  AOORESS 

University  of  Georgia 

Department  of  Chemistry 
Athens,  GA  30602 


It.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Office  of  Naval  Research 
Department  of  the  Navy 
Arlington,  VA  2'_i7 


14.  MONITORING  AGENCY  NAME  ft  ADORESV*'  dUltront  I 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


*.  RECIPIENT'S  CATALOG  NUMBER 


S.  TYRE  OF  REPORT  A  PERIOO  COVERED 

Technical  Report 


PERFORMING  ORO.  REPORT  NUMBER 


ONTRACT  OR  GRANT  NUMBERfaj 


N00014-84-K-0365 


12.  REPORT  DATE 

H/5/86 


IS.  NUMBER  OF  PAGES' 
16 


Controlling  Offleo)  I  IS.  SECURITY  CLASS,  (at  title  report) 


tSa.  'CLASSIFICATION/ DOWNGRADING 

schedule 


IS.  DISTRIBUTION  STATEMENT  (ol  i hit  Roport) 


This  document  has  been  approved  for  public  release  and  sale;  its  distribution  is  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  It.  obttroet  tntorod  In  Block  30,  II  dllloront  tram  Roport) 


IS.  SUPPLEMENTARY  NOTES 


To  be  published  in. 


IS.  KEY  WOROS  (Continue  on  rororto  tldo  It  nocoeeory  mnd  Identity  by  block  number) 

spin  dynamics,  Langevin  Equation,  XY  model  and  Ising  model 


20.  ABSTRACT  (Continue  on  r»v«ra«  cldo  It  no  cm  ••  ary  Old  Identity  by  block  numbot) 

We  investigate  the  dynamics  of  the  one-dimensional  S=%  isotropic  XY  model 
and  transverse  Ising  model  in  the  high  temperature  limit  by  using  the  method  of 
recurrence  relations.  We  obtain  the  relaxation  functions  as  well  as  some  Brownian 
analogs  of  a  generalized  Langevin  equation  for  a  tagged  spin  S  ?  in  these  models, 
namely,  the  memory  functions  and  the  random  forces.  We' find  that  the  realized 
dynamical  Hilbert  spaces  of  the  two  models  have  the  same  structure,  which  leads 
to  similar  dynamical  behavior  apart  from  a  time  scale. 


DO  ,iSTw  U73 


EDITION  OF  1  NOV  «•  IS  OBSOLETE 
S/N  0102- LF- 014- 4401 


Unclassified _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dote  tntorod) 

86  11  20  046 


m\ 


OFFICE  OF  NAVAL  RESEARCH 


Control  N00014-84-K-0365 
TECHNICAL  REPORT  NO.  32 

RELAXATION  FUNCTIONS,  MEMORY  FUNCTIONS,  RANDOM  FORCES,  AND  ERGODICITY 
IN  THE  ONE-DIMENSIONAL  SPIN-*  XY  AND  TRANSVERSE  ISING  MODELS 


by 

Joao  Florencio  Jr.  and  M.  Howard  Lee 


Prepared  for  publication  in 
Phys.  Rev.  B 


University  of  Georgia 
Department  of  Physics 
Athens,  Georgia  30602 


November  5,  1986 

Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  United  States  Government. 

This  documer.t  has  been  approved  for  public  release 
and  sale;  its  distribution  is  unlimited. 


RELAXATION  FUNCTIONS,  MEMORY  FUNCTIONS,  RANDOM  FORCES,  AND  ERGODICITY 
IN  THE  ONE-DIMENSIONAL  SPIN-J5  XY  AND  TRANSVERSE  ISING  MODELS 


JOAO  FLORENCIO  JR.  and  M.  HOWARD  LEE 
DEPARTMENT  OF  PHYSICS,  UNIVERSITY  OF  GEORGIA,  ATHENS,  GEORGIA  30602 


0  »  “  r 


mics  of 


J 


the  one-dimensional  S=%  isotropic 


•  f  r  LsO ‘J  A  •?  rT" 

We  investigate^the  dynai 
XY  model  and  transverse  Ising  model  in  the  high  temperature  limit 
by  using  the  method  of  recurrence  relations.  -We  obtain  the  relaxation 
functions  as  well  as  some  Brownian  analogs  of  a  generalized  Langevin 


equation  for  a  tagged  spin'll  in  these  models,  namely,  the  memory 

'-r  if  A- 

functions  and  the  random  forces.  -We-find- that  the  realized  dynamical 
Hilbert  spaces  of  the  two  models  have  the  same  structure,  which 
leads  to  similar  dynamical  behavior  apart  from  a  time  scale.  Based  on 

,t/’r 

the  infinite  dimensionality  of  these  Hilbert  spaces  we- also  conclude/^ 
that  is  ergodic  in  both  models.  <^- - - 


I.  INTRODUCTION 


The  one-dimensional  spin-%  XY  model  has  been  of  considerable 
theoretical  interest  in  recent  years  as  a  solvable  many-body  system1. 
The  Hamiltonian  of  this  model  is  given  by 


H  =  2  I  (JXS*S*+1  +  /S^+1)  -  B  z  S*  , 


(1.1) 


where  s“  are  spin  operators,  Ja  are  the  coupling  constants,  and  B  is 

an  external  magnetic  field.  Periodic  boundary  conditions  are  imposed, 

so  that  $“+1  =  S“  ,  where  N  is  the  total  number  of  spins  and  a  «  x,y,or  z 

In  this  paper  we  are  concerned  with  two  particular  cases  of  this 

Hamiltonian,  namely,  the  isotropic  XY  model  (XY)  for  which  J^J^J,  B=0, 

and  the  transverse  Ising  model  (TI)  where  JX=J,  J^O,  and  B=J. 

Although  the  equilibrium  properties  of  these  systems  are 
1  2 

well  known  *  ,  their  dynamical  behavior  is  less  well  understood. 

There  are  exact  results  for  the  longitudinal  time-dependent  spin 

3 

correlation  functions  due  to  Niemeijer  ,  and  for  the  transverse 

correlation  functions  in  the  limit  of  high  temperature  obtained  by 
4  5 

Brandt  and  Jacoby  and  also  by  Cape!  and  Perk  .  The  transverse  spin 
correlation  functions  for  both  the  XY  and  TI  cases  in  the  high  temperature 
limit  are  found  to  be 


(1.2) 


where  a  =  J2  for  the  XY  case,  and  J2/2  for  the  TI  case.  -  The  underlying 
reason  why  the  two  cases  have  the  same  time  dependence  has  not  been 
recognized  so  far. 


In  this  paper  we  use  the  method  of  recurrence  relations  due  to 
Lee®  to  study  the  time  evolution  of  both  the  isotropic  XY  model  and 
the  transverse  Ising  model  in  the  high  temperature  limit.  This  method 
allows  one  to  obtain  a  detailed  yet  rigorous  description  of  the 

dynamics  in  such  systems.  It  has  been  applied  to  some  spin  models^, 

8  9 

to  an  electron  gas  ,  to  a  classical  harmonic  oscillator  chain  ,  and 

to  the  study  of  velocity  autocorrelation  functions10. 

In  the  method  of  recurrence  relations  the  time  evolution 
of  a  dynamical  variable,  e.g.  S^t),  is  described  as  an  orthogonal 
expansion  in  a  properly  defined  Hilbert  space,  where  the  time  dependency 
is  placed  on  the  expansions  coefficients.  By  inspecting  the  relative 
norms  of  the  dynamical  Hilbert  space  for  each  of  the  XY  and  TI  cases 
we  can  readily  see  why  these  systems  are  dynamically  equivalent  in 
the  high  temperature  limit.  With  relatively  little  effort  we  recover 
the  transverse  correlation  function  (1.2)  for  these  systems.  We  would 
also  like  to  point  out  that  these  systems  are  dynamically  equivalent 
to  the  spin  Van  der  Waals  model  studied  by  Lee,  Kim,  and  Dekeyser^. 

We  conjecture  that  there  may  also  be  other  systems  with  similar 
dynamical  behavior,  for  which  their  respective  dynamical  Hilbert 
spaces  have  the  same  geometry,  that  is,  the  same  dimensionality 
and  also  the  same  relative  norms  of  basis  vectors. 

To  obtain  further  insight  in  the  time  evolution  of  these  systems 
we  calculate  some  Brownian  analogs  of  a  generalized  Lanqevin  equation 
for  the  spin  variables,  namely,  the  spin  memory  function  and  the 
spin  random  force.  We  also  discuss  ergodicity  in  these  models 
based  on  the  dimensionality  of  the  realized  Hilbert  spaces  of  the 
dynamical  variables  of  interest. 


The  arrangement  of  this  paper  is  as  follows.  In  Sec.  II  we  review 
the  method  of  recurrence  relations  as  well  as  its  connections  to 
a  generalized  Langevin  equation.  In. Sec.  Ill  that  method  is  applied 
to  the  dynamical  behavior  of  the  isotropic  XY  model  and  the  transverse 
Ising  model.  Correlation  functions,  relaxation  functions,  memory  functions, 
and  random  forces  are  then  obtained.  Finally,  in  Sec.  IV  we  summarize 
our  results  and  discuss  ergodicity  in  these  systems. 
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II.  METHOD  OF  RECURRENCE  RELATIONS 

Consider  a  one-dimensional  N  spin-**  system  described  by  a 
Hamiltonian  H.  The  time  evolution  of  an  operator  G  is  given  formally  by 


G(t)  =  eltL  G(0) , 


where  L  is  the  Liouville  operator  for  the  system,  defined  by 


(2.1) 


ILf  =  [f,H]  =  f H  -  Hf. 


Equation  (2.1)  can  also  be  expressed  by  the  expansion' 


d-1 

G(t)  *  I  a  (t)f  , 


(2.2) 


(2.3) 


where  fy  are  basis  vectors  of  a  Hilbert  space  A  of  d  dimensions. 
The  positive  definite  scalar  product  in  A  is  defined  in  the  high 
temperature  limit  (T=»)  as 


(A,B)  =  ~  Jr  AB+ 


(2.4) 


where  Z  =  2U  is  the  partition  function  of  the  system  in  this  limit. 

By  choosing  fQ  s  G(0)  it  follows  that  the  remaining  basis  vectors 


f  can  be  generated  by  the  following  recurrence  relation  (RRI): 


f  ,,  =  i Lf  +  A  f  ,  ,  v>0  , 

v+1  v  v  v-1  — 


(2.5) 


m 


5 


where 


A 

v 


v>l  , 


(2.6) 


are  the  relative  norms  of  the  basis  vectors,  and  by  definition 
f_J  =  0,  Aq  =  1. 

The  coefficients  av(t),  which  are  also  the  relaxation  functions, 
satisfy  a  second  recurrence  relation  (RRII): 


Av+iav+1(t)  3  -ay(t)  +  av-1(t)  ,  \>>0  ,  (2.7) 

.  da  (t) 

where  a  (t)  =  — —  ,  a  .  =  0.  Notice  that  due  to  the  initial  choice 
v'  '  dt  -1 

fQ  *  G(0) ,  it  follows  from  Eq.(2.3)  that  8q(0)  =  1,  and  av(0)  =  0 
for  y>l.  The  complete  time  evolution  of  G(t)  can  thus  be  determined 
by  using  RRI  and  RRII. 

A  generalized  Langevin  equation  for  the  operator  G(t),  which  is 
formally  equivalent  to  the  Heisenberg  equation  of  motion,  is  given  by 


+  /  df*(t-t')G(f)  =  F(t)  ,  (2.8) 

ac  o 

where  <j>  is  the  memory  function  and  F  the  random  force.  Both  <{>  and  F 
can  be  readily  obtained  as  follows.  The  random  force  is  given  by 


(2.10) 


t 

a  (t)  =  /  dt'b  (t-t')aQ(t'),  v>l. 

v  0 

The  memory  function  is  given  simply  by  <j>(t)  =  b^( t) -  The  remaining  b^'s 
that  is,  b2,  b3,  ...  ,  are  the  2nd  memory  function,  3rd  memory  function, 
etc.  The  reader  is  referred  to  the  original  formulation  of  the  method 
of  recurrence  relations  for  the  detailed  derivations  of  the  relations 
contained  in  this  section6. 
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III.  DYNAMICS  OF  THE  XY  AND  TI  MODELS 


In  this  section  we  apply  the  method  of  recurrence  relations 


described  in  Sec.  II  to  investigate  the  dynamical  properties  of 


the  XY  and  TI  models  at  T  s  ».  Consider  first  the  XY  case,  taking 


S^  as  the  dynamical  variable  of  interest.  The  time  evolution  of  S^(t) 

J  J 


is  given,  according  to  Eq.(2.3),  as 


S*(t)  =  z  a  (t)f  , 


(3.1) 


where  f0  =  S^O)  =  Sx  By  using  RRI  we  obtain 
fl  '  2J<Sj-lSjZ  +  SjSj+l>- 

f  =  -  4j2(sx  s2  S2  -  SX  SySy  +  2Sy  SxSy  -  Sy  SySx  +  S2S2  Sx  1 
f2  1  j-2bj-rj  bj-i  jbj+i  j-ijj+i  j-i  j  j+i  j  j+i  j+2; 

f?  =  2J3(-4Sy  -S2  ?S2  .S2  +  4Sy  ?S2  ,sysy  -  8SX  ?SZ  ,SySX  . 

3  '  j-3  j-2  j-1  j  j-2  j-1  j  j+1  J-Z  J-l  j  j+1 


+  12SX  9S2  ,SxSy  .  +  3S^  ,SZ  .  +  3S2  ,Sy  .  +  12Sy.  .S^S2  Sx 
J-2  j-1  j  j+1  j-1  j+1  j-1  j+1  j-1  j  j+1  j+2 


-  8SX  SySZ  SX  +  4Sy  SyS2  Sy  -  4S2S2  S2  Sy  1 

abj-ibjbj+ibj+2  +  4bj-rj5j+rj+2  4bjbj+r  j+2bj+3J » 


f  =  8J4ISX  SX  SX  +  SX  SX  Sx  -  2(Sy  SX  Sy  +  ^y  SX  Sy 

r4  00  (  j-2  j-rj+i  bj-i  j+ibj+2  Abj-2bj-ibj+i  -j-r j+rj+2 


-  sx  s2  s2  s2  s2  +  sx  sz  sz  sysy  +  sy  sysz  s2  sx 
bj-4bj-3bj-2 bj-ibj  bj-3bj-2bj-rjbj+i  j-i  j  j+i  j+2bj+3 


S2SZ  Sz  SZ  SX  )  +  3(S*  s*  +  s*  s*  sA  -  s4  sAs4  ) 

j  j+i  j+2^j+3bj+4;  +  J^j-2bj-rj+i  bj-r j+rj+2  bj-i  jbj+r 


;x  <.y  pX 


y  ,y  pX 


-Z  CXCZ 


"  4(Sj-2Sj-lSj+l  +  Sj-lSj+lSj+2)  "  6tSj-3Sj-2Sj-lSjbj+l 

+  sx  <>z  sysz  ys  +  Sy  <;z  cy qz  sx  +  sx  sysz  s2  sy  ) 
bj-2bj-lbjbj+lbj+2  bj-2bj-lbjbj+lbj+2  j-1  j  j+1  j+2bj+3; 


■"'••Nr 


I 


♦  8^Sj-3Sj-2Sj-lSjSj+l  +  Sj-lSjSj+lSj+2Sj+3^ 


-  12SX  S2  SXS2  SX  1 
12bj-2^j-r  r j+i  j+2'  * 


(3.2) 


Thus  an  excitation  of  at  t=0  will  propagate  through  the  chain 

according  to  Eqs.(3.1  and  2).  The  basis  vectors  f  correspond  to 

excitations  of  clusters  of  spin.  There  is  a  spatial  hierarchy  in 

the  f  's  in  the  sense  that  as  \>  increases,  so  do  the  boundaries  of 
v 

the  region  in  the  chain  within  which  clusters  of  spins  are  excited. 

As  we  shall  see  later,  there  is  also  a  hierarchy  in  the  time  sequence 
in  which  the  vectors  f  are  excited,  so  that  a  true  propagation  of 
the  initial  excitation  throughout  the  chain  does  indeed  take  place. 

In  addition,  the  length  of  the  clusters  generally  increases  with  v. 

Notice  also  the  appearance  of  "disconnected"  clusters  in  vectors  f^ 
of  higher  dimensions  of  A  consisting  of  groups  of  spins  in  which 
at  least  one  of  the  spins  is  separated  from  the  others,  e.g.  Sj_2Sj_i^j+i 
in  f4. 

The  relative  norms  a  are  easily  obtained  from  Eqs. (2.4,6)  and 

V 

(3.2).  We  find  A^  =  A,  a2  =  2a,  a^  =  3a,  etc.,  where  A  =  2J2. 

u  7 

The  quantity  A2  is  referred  to  in  the  literature  as  the  basal  frequency  . 
In  the  thermodynamic  limit  there  is  no  upper  bound  for  v,  so  that 
the  dynamical  Hilbert  space  A  has  infinite  dimensions  (d  =  <*>) , 
and  the  relative  norms  are  given  by 


A  =  vA , 

V 


x,  2,  3,  ... 


(3.3} 


A  relation  similar  to  this  was  also  obtained  by  Lee,  Kim,  and  Dekeyser 
in  their  study  of  the  dynamics  of  the  spin  Van  der  Waals  model.'7 


The  recurrence  relation  (2.8)  now  reads 


(v+l)aav+1(t)  -  -av(t)  +  a^_^(t),  v  0*  1,  2,  ...  (3.4) 


It  is  satisfied  by 

a  (t)  =  ~t  exp(-%At2)  . 

V  V  • 


(3.5) 


These  are  the  so-called  relaxation  functions. 

.L 

In  terms  of  normalized  basis  vectors  F  =  (f  ,f  )  f  , 

V  V  V  V 

the  time  evolution  of  S^(t)  is  given  by 

J 


v=0 


A  (t)F 

v  V 


(3.6) 


where 


A(t)  =  .  (3.7) 

V  2  (v'P 

This  quantity  satisfies  the  Bessel  equality 

I  A2(t)  =  I  .  (3.8) 

— n  ^  ** 

v-0 

That  is,  the  length  of  the  vector  S*(t)  in  the  dynamical  Hilbert  space 

0 

is  an  invariant  of  time.^’^  In  Fig.l  we  show  the  time  dependent 
probabilities  A2(t)  (normalized  to  k)  in  early  stages  where  S*(t)  samp] 
the  space  of  the  lower  basis  vectors.  Notice  that  due  to  the 
pre-exponential  factor  in  A  (t),  the  basis  vectors  corresponding  to 
the  lower  dimensions  of  are  more  likely  to  be  initially  excited 
before  those  of  higher  dimensions.  In  addition,  each  of  the 


probabilities  A2(t)  decay  in  a  Gaussian  manner,  so  that  a  somewhat 
localized,  yet  spreading,  group  of  basis  vectors  are  being  excited 
as  time  progresses.  This  corresponds  to  a  propagation  of  the  initial 
excitation  throughout  the  spin  system.  This  propagation,  however, 
is  to  be  understood  In  a  quantum  mechanical  context,  since  the  Bessel 
equality  (3.8)  holds  for  all  times.  That  is  because  actually  all  the  f 
are  excited  for  t>0  due  to  the  instantaneous  nature  of  the  coupling 
constant.  However,  at  any  given  t  the  probability  of  excitation  Is 
significant  only  for  a  fraction  of  the  f^'s.  On  the  other  hand,  it  can 
be  seen  from  Eq -(3.1)  that  the  maxima  of  the  amplitudes  A2(t)  are 
located  at  v  -  At2.  In  the  spin  chain,  that  corresponds  to  the 
boundaries  j±v  of  the  region  in  which  the  significant  excitations 
are  taking  place. 

The  transverse  correlation  function  can  now  be  obtained  as  follows 


<S*(t)S*(0)>  =  (S*(t),S*)  =  £a0(t) 

=  i  exp(-*sAt2) 

=  ^  exp(-02t2)  ,  (3.9) 


4  5 

which  agrees  with  the  results  found  in  the  literature.  * 

Consider  now  the  TI  case.  The  time  evolution  of  the  operator 
is  represented  by 

J 


11 


V  V 

where  f'0  =  S ^ (0)  =  ,  and  primed  variables  are  used  to  distinguish 

the  notation  here  from  that  of  the  XY  case.  By  using  RRI  we  obtain 


f',  - 


f'o  = 


f'o  = 


f.  = 


Jsj  • 


-2J2<sj-isj  +  sjsj+i>  • 

-2J3(sj-isj  +  4sj-isjsj+i  +  sjsj+i*  * 
4J**(SX  ?SZ  .S?  -  2S5}  ,S^SX  ,  -  3SX  .$ySy 

J-H  J-l  J  J-l  J  J+l  j-l  j  j+l 

ZSZ  SX  ) 
j  j+l  j+2J 


3Sy  SySx 
j-l  j  j+l 


+  S.;S  •  ..S . .«)  ,  etc. 


(3.11) 


In  spite  of  the  basis  vectors  f'  here  are  different  from  the  basis 

vectors  f  of  the  XY  case,  the  relative  norms  assume  the  same  form, 
v 

namely, 

A'v  =  va'  ,  v  -  1,  2,  3 .  (3.12) 

where  a'v  =  J2.  This  leads  to  the  same  dynamical  behavior  as  with  the 
XY  model  and  the  spin  Van  der  Waals  model.  In  particular,  the  spin 
correlation  function  for  the  TI  case  is  found  to  be 


<Sj(t)Sj(0)>  =  ^  exp(-JsJ2t2)  ,  (3.13) 

4  5 

which  also  agrees  with  previous  results.  *  Notice  that  the  basal 
frequencies  a*5  and  a'1*  of  these  models  can  be  made  equal  by  rescaling 

the  coupling  constants.  Conversely,  one  can  say  that  the  XY  model 

I 

and  the  TI  model  are  dynamically  equivalent  with  respect  to  the  time 
evolution  of  S^(t). 

’  J 

i 

i 

I 


i 


Consider  now  the  Brownian  analogs  of  the  generalized  Langevin 
equation,  Eq. (2.8) .  The  spin  random  forces  for  the  XY  and  TI  cases 
are  given  respectively  by 

F(t)  »  t  b  (t)f  (3.14a) 

V=1 

and 

F'(t)  =  Z  b'v(t)f'v  ,  (3.14b) 

Vs*  1 

with  f  and  f'  given  by  (3.2)  and  (3.11).  The  memory  functions  b^(t) 
can  be  expressed  as 


b  (t)  -  z  al5(m-''+1)  Cl  t1"  .  (3.15) 

v  m=0  * 

with  a  similar  expression  for  b'^(t)  involving  a*.  The  coefficients 
are  obtained  recursively  by 


2p 

E 

r=0 


2PpI 


1 

2n+v-r 


C 


v 

2(n-p)+v-l 


f-ln 

-i—1 —  ,  n=0,l,2,... 
2nn!vl 


(3.16) 


The  memory  function  $(t)  of  the  XY  model  is  found  to  be 


♦(t)  =  b2(t)  =  1 


2At2  .  10A2t4  74A3t6  .  706A4t8  8612A5t10  , 

T T~  “43  S3  83~  TS3 


(3.17) 
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A  similar  expression  for  the  memory  function  of  the  TI  model  is 
obtained  from  Eq. (3.17)  by  the  replacement  A  ■+•  a'.  The  remaining 
memory  functions  bv(t)  were  first  calculated  by  Lee  et  al^,  and 
their  results  are  reproduced  here  in  Fig. 2.  In  that  figure,  A  is 
taken  to  be  unity  for  convenience.  From  the  generalized  Langevin 
equations  one  can  see  now  that  S^(t)  evolves  in  time  modulated  by 
the  memory  function  <j>.  The  random  force  F,  with  components  along 
the  higher  dimensions  of  the  dynamical  Hilbert  space  (v  s  1  -*■  «) 
acts  to  pull  the  time  evolution  away  from  the  basal  plane  onto 
higher  reaches  of>A  . 


IV.  CONCLUSIONS 


We  investigated  the  time  evolution  of  a  tagged  spin  variable  S*(t) 

J 

in  the  one-dimensional  Sa%  isotropic  XY  model  and  transverse  Ising  model 
at  infinite  temperature  by  using  the  method  of  recurrence  relations. 

This  method  lends  itself  to  a  detailed  description  of  the  dynamical 
behavior  of  these  systems,  whether  by  looking  at  the  propagation  of 
the  excitations  along  the  spin  chains,  or  by  examining  the  effects 
of  the  memory  funcions  and  the  random  forces  of  a  generalized  Langevin 
equation.  The  time  evolution  of  S.(t)  is  given  as  an  orthogonal 

J 

expansion  in  a  properly  defined  Hilbert  space  in  each  of  these  systems. 

The  relative  norms  of  the  basis  vectors  in  each  of  the  XY  and  TI  models 
have  the  same  structure,  resulting  in  similar  dynamical  behavior  for  S^t) 
in  both  cases.  We  obtained  expressions  for  the  transverse  spin  correlation 
functions,  relaxation  functions,  memory  functions,  and  random  forces 
for  these  systems.  These  results  are  valid  in  the  thermodynamic  limit 
(N  *  •). 

For  a  finite  system,  there  is  always  an  upper  bound  for  the 
dimensionality  d  of  the  realized  Hilbert  space  and,  as  a  consequence, 
only  periodic  solutions  are  admissible.  In  that  case  the  system  will 
eventually  return  to  its  initial  state  with  an  excitation  localized, 
say,  at  site  j.  Thus,  for  a  finite  system,  there  can  never  be 
an  equi partition  of  energy  and  the  system  is  non-ergodic.  On  the 
other  hand,  for  infinite  systems  the  dimensionality  of  the  Hilbert 
space  has  no  a  priori  constraints  on  its  upper  bound.  The  dimensionality 
of  the  Hilbert  space  depends  on  the  model  as  well  as  on  the  dynamical 
variable  under  consideration.1* 
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As  we  have  seen,  the  Hilbert  space  of  S.(t)  has  infinite  dimensions 

J 

in  both  XY  and  TI  cases  discussed  here.  This  leads  to  a  propagation 

of  an  excitation  throughout  the  spin  chain,  whose  propability 

amplitudes  decay  with  time  in  a  Gaussian  fashion.  The  system 

under  consideration  never  returns  to  its  initial  state,  thus  the 

process  is  irreversible.  As  t  «,  having  an  infinite  number 

of  basis  vectors,  ( t)  must  sample  all  reaches  of  the  Hilbert  space 

and  it  is  ergodic  in  both  the  XY  and  TI  models.  This  conclusion  is  in 

12 

accordance  with  the  results  of  Perk  et  al  for  the  XY  case.  In  the 

case  where  the  dynamical  variable  of  interest  is  a  conserved  quantity, 

the  dimensionality  of  its  Hilbert  space  is  d*l,  and  that  variable 

is  non-ergodic.  That  is  the  case  with  the  magnetization  in  the 
12  13 

XY  model  *  .  We  are  currently  investigating  the  formal  aspects 

of  the  connection  between  ergodicity  and  dimensionality  of  realized 
dynamical  Hilbert  spaces.  Our  results  will  be  reported  in  a  forthcoming 
paper. 
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FIGURE  CAPTIONS 


Fig.  1  Probabilities  A2(t)  vs  time»S^(t)  samples  the  basis  vectors  f  . 

k 

The  time  is  given  in  units  of  the  basal  frequency  er. 

Fig.  2  Memory  functions  bv(t)  vs  time  for  both  XY  and  TI  models. 

The  basal  frequencies  A5  and  are  taken  to  be  unity  in 
both  axes. 


